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We present the results of Monte Carlo simulation studies of the capability of the proposed NEMO telescope to 
detect TeV muon neutrinos from Galactic microquasars. In particular we determined the number of the detectable 
events from each known microquasar together with the expected atmospheric neutrino and muon background 
events. We also discuss the detector sensitivity to neutrino fluxes expected from microquasars, optimizing the 
event selection in order to reject the atmospheric background, and we show the number of events surviving the 
event selection. 



1. INTRODUCTION 

Microquasars are Galactic X-ray binary sys- 
tems which exhibit relativistic jets, observed in 
the radio band pQ. Several authors propose mi- 
croquasar jets as sites of acceleration of charged 
particles up to energies of about 10 16 eV, and of 
high energy neutrino production. According to 
present models, neutrinos could be produced both 
in p7 (213) and pp interaction scenarions (4151617) . 

The aim of this paper is to study the possibil- 
ity to detect neutrinos from known microquasars 
with the proposed NEMO-km 3 telescope In 
particular, for each microquasar we simulated the 
number of expected events according to present 
models. We also simulated the expected back- 
ground from atmospheric neutrinos and muons. 

We applied an event selection in order to reject 
the atmospheric background. The event selection 
is a combination of different criteria, chosen opti- 
mizing the detector sensitivity to neutrino fluxes 
coming from each microquasar. Eventually we 
calculated the expected number of events surviv- 
ing the selection referring to 1 year of data taking. 

2. DETECTOR LAY-OUT 

The NEMO-km 3 telescope, simulated in this 
paper, is a square array of 9 x 9 towers with a dis- 
tance between towers of 140 m. In this configura- 
tion each tower hosts 72 PMTs (with a diameter 
of 10"), namely 5832 PMTs for the whole detec- 
tor with a total geometrical volume of ~ 0.9 km 3 . 



We considered an 18 storey tower; each storey is 
made of a 20 m long beam structure hosting two 
optical modules (one downlooking and one look- 
ing horizontally) at each end (4 OMs per storey). 
The vertical distance between storeys is 40 m. A 
spacing of 150 m is added at the base of the tower, 
between the anchor and the lowermost storey. 

The detector response is simulated using the 
simulation codes developed by the ANTARES 
Collaboration 9 . In the simulation codes, the 
light absorption length, measured in the site of 
Capo Passero (L a w 68 m at 440 nm ) , is 
taken into account. Once the sample of PMT 
hits is generated, spurious PMT hits, due to the 
underwater optical noise ( K decay), are intro- 
duced, with a rate of 30 kHz for 10" PMTs, corre- 
sponding to the average value measured in Capo 
Passero site. 

3. EXPECTED NEUTRINO EVENTS 
FROM MICROQUASARS 

For each microquasar, a number N to tai = 10 9 of 
interacting neutrinos was simulated in the energy 
range 1 - 100 TeV in order to reproduce the neu- 
trino fluxes calculated by Distefano et al. |3] • The 
generation spectral index is chosen to be X = 1, 
in order to guarantee a good event statistics at 
the highest energies. The simulated events are 
then weighted to the theoretical fluxes and the 
number of expected events is computed and re- 
ported in Tab. ^ Results refer to an integration 
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time At equal to the duration of the considered 
burst for the transient sources and to one year for 
the steady sources. 

For microquasar LS 5039, N tota i = 10 8 events 
were generated in the energy range 0.1-1 TeV to 
simulate neutrinos events according to the model 
proposed by Aharonian et al. 5 . In Tab. [21 
the number of muon events from LS 5039 is re- 
ported. The simulations refer to a power-law neu- 
trino spectrum, dn v /de v oc e~ r , with energy cut- 
off e™ ax = 10 TeV and 100 TeV, and T = 1.5 and 
2, respectively. For the four combinations of pa- 
rameters r and £™ ax , an average neutrino energy 
flux fl h = 10- w erg/cm 2 s {e v > 0.1 TeV) was 
simulated. 

For microquasar Cygnus X-3, N to tal — 10 8 
events were generated in the energy range 0.1 
- 1 TeV and N to tai = 10 9 events with energies 
100 - 1000 TeV to simulate Bednarek's theoreti- 
cal fluxes @|. In Tab. [3] we give the number of 
muon events from Cygnus X-3, considering dif- 
ferent values of the the spectral index k and for 
the Lorentz factor cut-off j^ ax of the spectrum of 
nuclei generating the neutrinos. The case of neu- 
trons injected by mono-energetic nuclei having a 
Lorentz factor ^<™ lm « 10 5 is also considered (see 
H] for details). 

Other models, present in literature, are not 
considered here because they refer to micro- 
quasars outside the field of view of the NEMO 
telescope |Hj or because the predicted neutrino 
fluxes are too low to be detected by a km 3 tele- 
scope 0. 

3.1. SIMULATION OF THE BACK- 
GROUND 

A sample of 7 • 10 9 atmospheric neutrinos have 
been generated using the ANTARES event gener- 
ation code, based on a weighted generation tech- 
nique 9 . The events were generated in the en- 
ergy range 10 2 -j- 10 8 GeV, with a spectral index 
X = 2 and a Air isotropic angular distribution. 
The events were then weighted to to the sum of 
the Bartol flux |llj and of prompt neutrino rqpm 
model 52 nux - So doing, we compute a number 
pa 4 • 10 4 of detected atmospheric neutrino events 
per year of data acquisition. 

Atmospheric muons are generated at the de- 



tector, applying a weighted generation technique. 
We generated a sample of N tota i = 2.5- 10 7 muons, 
in the energy range 1 TeV -j- 1 PeV, with a gen- 
eration spectral index X = 3. We also generated 
Ntotai = 4 • 10 7 events in the range 100 GeV - 
1 TeV, with a generation spectral index X = 1. 
Muons are generated with an isotropic angular 
distribution. The events are weighted to the 
Okada parameterization |13|. taking into account 
the depth of the Capo Passero site (D — 3500 m) 
and the flux variation inside the detector sensi- 
tive height (h « 900 m) . According to the Okada 
parameterization, the expected number of recon- 
structed muon events is about 4 • 10 8 per year. 
Our statistics cover only a few days. Considering 
that reconstructed events have a flat distribution 
in Right Ascension (RA), we can project the sim- 
ulated events in a few degrees bin ARA, centered 
in the source position. So doing, we get statistics 
of atmospheric muons corresponding to a time ^1 
year at all source declinations. 

At the reconstruction level, the fraction of 
background events is too high and exceeds the ex- 
pected signal also in the case of the most intense 
microquasars. As an example we present the 
case of microquasars SS443 ((5=+04° 58' 58" .0, 
RA=19 ;i ll m 49 s . 6), considering an observation 
time of 1 year and an angular bin of 3° x 3° around 
the source position in the sky. The result for 
SS433 is shown in Fig. [J]in which are represented 
Nfj, — 71 events from the source and = 2113 
from atmospheric background (only up-ward go- 
ing reconstructed tracks are considered). Source 
and background events are selected with a hit-or- 
miss procedure. This procedure consists in com- 
paring, for each event, the weight W even t with a 
number R randomly generated in the range [0, 1]. 
If R < W event the event is accepted, if not the 
event is rejected. 



4. SENSITIVITY TO MUON NEUTRI- 
NOS FROM MICROQUASARS 

4.1. Calculation of the detector sensitivity 

The detector sensitivity was calculated accord- 
ing to the Feldman and Cousins approach |14|. 
The 90% c.l. sensitivity to a neutrino flux com- 
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Table 1 

Expected number of neutrino induced muons from microquasars expected from the theoretical neutrino 
energy flux quoted by Distefano et al. |2| : Nf t xp is the total number of reconstructed muons from each 
microquasar during the time interval At, N™ is the expected number of selected muon neutrino events 
from the sources and N b is the atmospheric background events surviving the event selection and expected 
in 1 year of data taking. 



Steady Microquasars 


At (days) 


ivrexp 

V 


N m 




LS 5039 


365 


0.2 


0.1 


0.1 


Scorpius X-l 


365 


0.8 


0.2 


0.1 


SS433 


365 


205.0 


76.0 


0.1 


GX 339-4 


365 


185.0 


68.0 


0.1 


Cygnus X-l 


365 


2.0 


0.5 


0.1 


Bursting Microquasars 


XTE J1748-288 


20 


2.2 


0.8 


0.3 


Cygnus X-3 


3 


3.4 


0.8 


0.1 


GRO J1655-40 


6 


1.7 


0.6 


0.1 


GRS 1915+105 


6 


0.4 


0.1 


< 0.1 


Circinus X-l 


4 


0.2 


0.1 


0.1 


XTE J1550-564 


5 


< 0.1 


< 0.1 


< 0.1 


V4641 Sgr 


0.3 


< 0.1 -e- 4.6 


< 0.1-7-1.4 


0.1 


GS 1354-64 


2.8 


< 0.1 


< 0.1 


0.1 


GRO J0422+32 


l-r-20 


O.l-r-1.5 


< 0.1-7-0.4 


0.1 


XTE J1118+480 


30-7-150 


2.4-7-12.0 


1.0-7-4.8 


0.2 



Table 2 

Expected number N^ xp of neutrino induced muons from microquasar LS 5039 according to Aharonian et 
al. model 0, during an time interval of 1 year. We also give the number N™ from the source compared 
to the atmospheric background events N b surviving the event selection. 



r 


£ rnax (IfeV) 


N exp 




N b 

M 


1.5 


10 


7.0 


1.7 


0.2 


1.5 


100 


15.0 


4.9 


0.1 


2.0 


10 


4.3 


1.0 


0.3 


2.0 


100 


9.0 


2.6 


0.1 
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Table 3 



Expected number of neutrino induced events from microquasar Cygnus X-3 (Bednarek model 0]): N^ xp 
is the total number of reconstructed muons during an time interval of 1 year; iV™ and are the number 
of source and background events surviving the selection. 
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REG. II 
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V 


ivrexp 
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V 


2.0, 10 6 


2.0 


1.8 


0.4 


0.4 


0.1 


2.0, 10 7 


5.2 


5.0 


1.1 


1.1 


0.1 


2.5, 10 7 


0.3 


0.2 


< 0.1 


0.1 


0.1 


mono-energetic 


< 0.1 


4.5 


< 0.1 


0.7 


0.1 



source only source-!- background 




RA (h) RA (h) 

Figure 1. Maps of the reconstructed events from SS433 (left panel) during an observation time of 1 year. 
In the right panel background events are plotted together with the source events. 



ing from a point-like source is given by 
M 90 ( & ) 



fv, 



90 



N„ 



-fv,0, 



(1) 



where JI 90 (b) is the 90% c.l. average upper limit 
for an expected background with known mean 
value b and no true signal f v< o is an arbitrary 
point source flux inducing a mean signal iV M . 

4.2. Selection criteria for background re- 
jection 

The used reconstruction algorithm is a robust 
track fitting procedure based on a maximization 
likelihood method 9 . In this work, we used, as 



a goodness of fit criterion, the variable: 



A 



log(£) 



N 



DOF 



■0.1(N, 



camp 



i), 



(2) 



where log(£) /Ndof is the log-likelihood per de- 
gree of freedom {Ndof) and N comp is the total 
number of compatible solutions found by the re- 
construction program. In particular, events are 
selected if the variable A is greater than a given 
value A cut . This quality cut is here applied to- 
gether with other selection criteria as listed in 
the following: 

• the number of hits Nfa , used to reconstruct 
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the muon track, must be greater than a 
given value 

• the muon must be reconstructed with 
"in order to reject down-going 



,Qrec ^ „amax 

events; 



• only events reconstructed in a circular sky 
region centered in the source position and 
having a radius of Tun are considered. 

The optimal values of A cut , Nffi , $™ ax and r bm 
are chosen to optimize the detector sensitivity. 

4.3. Results 

The detector sensitivity was calculated for a 
livetime of 1 year, simulating a neutrino flux with 
spectral index T = 2 in the energy range 1-100 
TeV and taking into account both atmospheric 
neutrino and muon background. The study was 
carried out for each microquasar, since the sen- 
sitivity is a function of the source astronomical 
declination. 

In Tab. 0]we give the sensitivity for each mi- 
croquasar, and the corresponding i?™" 1 , A cut and 
run- The dependence of the sensitivity on the 
source declination is shown in Fig. [3 

In Tab. |S]are summarized the detector sensitiv- 
ities for microquasar LS 5039, assuming neutrino 
fluxes with spectral indices T — 1.5 and 2, in the 
energy range 0.1 TeV and e" lax — 10 and 100 
TeV. 



5. EXPECTED NUMBER OF EVENTS 

We applied the event selection cuts given in 
Tab. 0]and[Slto determine the expected number of 
events from theoretical models considered in this 
paper. In Tab. ^we give the number of selected 
neutrino events from each microquasar, according 
to the neutrino fluxes given by Distefano et al. [H]. 
In the same table we compare these results with 
the atmospheric neutrino and muon background 
in 1 year of data taking. In this analysis, it is as- 
sumed that transient sources cause one burst per 
year, i.e. the number of source events produced 
in the interval At is related to 1 year observation 
time. 

The search for neutrino events in coincidence 
with microquasar radio outbursts could restrict 




<5 (deg) 

Figure 2. NEMO-km 3 sensitivity to neutrinos 
from microquasars versus source declination, for 
a livetime of 1 year. The worsening of the sen- 
sitivity with increasing declination is due to the 
decrease of the source visibility. 



the analysis period to the flare duration At, re- 
ducing the background. Such an analysis tech- 
nique, already used by AMANDA [T3], can im- 
prove the detector sensitivity to neutrinos from 
transient sources. Referring to the bursts con- 
sidered in Tab. ^ and integrating over the time 
interval At, we expect about 10 -3 background 
events per burst. 

For the microquasar LS 5039 we considered the 
flux predicted by Aharonian et al. The ex- 
pected number of selected events is given in Tab. 
121 the comparison with the atmospheric back- 
ground shows that an evidence could be expected 
in a few years of data taking. 

In Tab. |31 we quote the expected number of 
events from microquasar Cygnus X-3, according 
to the Bednarek model 0]. Since the model pre- 
dicts a neutrino flux with a spectral index close to 
2 and since the main event contribution is in the 
energy range 1-100 TeV, we used the same event 
selection parameters quoted in Tab. 0] Except 
the case of re = 2.5 and 7™ 1 = 10 7 , the num- 
ber of events exceeds the expected background 
and the signal could be detected in a few years of 
data taking. 
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Table 4 



Detector sensitivity to neutrinos from microquasars: the sensitivity f v 
spectrum in the energy range 1-100 TeV, for a detector live time of 
of {)™ ax , A cut and r^n and the source declination 5 are also given. 


is calculated for an e v 2 neutrino 
1 year. The corresponding values 


Steady Microquasars 


$rnax ( deg ) 


A CM t 


run (deg) 


f v (erg/cm 2 s) 


LS 5039 


101 


-7.3 


0.9 


6.5- 10" 11 


Scorpius X-l 


104 


-7.7 


0.7 


5.8- 10~ n 


SS433 


115 


-8.0 


0.8 


5.7- 10- 11 


GX 339-4 


96 


-7.4 


0.5 


4.7- 10~ n 


Cygnus X-l 


103 


-7.5 


0.7 


9.0- 10- 11 


Bursting Microquasars 


XTE J1748-288 


102 


-7.6 


0.9 


5.4- 10~ n 


Cygnus X-3 


101 


-7.3 


0.8 


1.1 • io- 10 


GRO J 1655-40 


101 


-7.4 


0.7 


5.2- 10~ n 


GRS 1915+105 


100 


-7.4 


0.8 


7.4- 10~ n 


Circinus X-l 


90 


-7.3 


0.9 


4.2- 10~ n 


XTE J1550-564 


90 


-7.1 


0.9 


4.4- 10~ n 


V4641 Sgr 


102 


-7.4 


0.9 


5.6- IO" 11 


GS 1354-64 


90 


-7.5 


1.0 


3.8- IO" 11 


GRO J0422+32 


103 


-7.5 


0.8 


8.7- IO" 11 


XTE J1118+480 


102 


-7.5 


0.7 


1.1 • io- 10 



Table 5 

Detector sensitivity flux to neutrinos from LS 5039: the sensitivity /„ is calculated for a e~ r neutrino 
spectrum in the energy range from 0.1 TeV up to e™ ax and for a detector live time of 1 year. The 
corresponding values of $™ ax , A CMt and run are also given. 



£ C (TeV) $™ ax (deg) Acut r bm (deg) /„ (erg/cm 2 s) 

1.5 10 101 -7.3 1.0 1.5 -IO" 10 

1.5 100 101 -7.3 0.9 5.2 -10- 11 

2.0 10 97 -7.4 1.0 2.7 -IO" 10 

2.0 100 101 -7.3 0.9 9.6 -IO" 11 
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6. CONCLUSIONS 

The possibility to detect TeV neutrinos from 
Galactic microquasars with the proposed NEMO- 
km 3 underwater Cerenkov neutrino telescope has 
been investigated. 

A Monte Carlo was carried out to simulate the 
expected neutrino-induced muon fluxes produced 
by point-like sources and by atmospheric neu- 
trinos. The expected atmospheric muon back- 
ground was also simulated. 

We computed the detector sensitivity for each 
microquasar, optimizing the event selection in or- 
der to reject the background. Finally, we applied 
the event selection and calculated the number of 
surviving events. 

Our results show that, assuming reasonable 
scenarios for TeV neutrino production, the pro- 
posed NEMO telescope could identify micro- 
quasars in a few years of data taking, with a dis- 
covery potential for at least few cases above the 
5cr level, or strongly constrain the neutrino pro- 
duction models and the source parameters. 
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